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The vortex density of a rotating superfluid, divided by its particle mass, dictates the superfluid's 
angular velocity through the Feynman relation. To find how the Feynman relation applies to su- 
perfluid mixtures, we investigate a rotating two-component Bose-Einstein condensate, composed of 
bosons with different masses. We find that in the case of sufficiently strong interspecies attraction, 
the vortex lattices of the two condensates lock and rotate at the drive frequency, while the superfiu- 
ids themselves rotate at two different velocities, whose ratio is the ratio between the particle mass 
of the two species. In this paper, we characterize the vortex-locked state, establish its regime of 
stability, and find that it surives within a disk smaller than a critical radius, beyond which vortices 
become unbound, and the two Bose-gas rings rotate together at the frequency of the external drive. 

PACS numbers: 



After the first experimental realization of Bose- 
Einstein condensates (BECs) of alkali atoms, their study 
has experienced enormous advancements [l| . Among the 
major threads of investigation in BECs has been the 
study of vortices both experimentally 0, [^, 0] and theo- 
retically . It is known from the classic works of Onsager 
and Feynman [^, that superfluids rotate by nucleating 
vortices. When there are several vortices present, they 
form a triangular Abrikosov vortex lattice [8| , with den- 
sity given by 
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where m is the mass of a constituent boson and ft is the 
rate at which the superfluid - which rotates with the vor- 
tex lattice - is being rotated (see, for instance, Ref. [§]). 
The so-called "Feynman relation" ([1]) states that, on av- 
erage, a uniform superfluid rotates like a rigid body. It 
has been shown that corrections to Eq. ([T]) due to the 
typical experimental situation of nonuniform superfluid 
density resulting from a harmonic trap are small [lo| 
but experimentally observable [llj. Vortex physics be- 
comes much more intriguing in multi-component BECs. 
So far, the investigation of vortex lattices in multi- 
component BECs utilized different hyperfine levels of 
the constituent atoms to obtain multi-component con- 
densates (e.g., 12, [l^). Thus the mass of all conden- 



sate components was identical, and the generalization of 
Eq. ^ straightforward. 

In this paper we ask what are the consequences of the 
Feynman condition, Eq. ([T]), in a system of interact- 
ing two-component rotating condensates with different 
masses. We find that for sufficiently large attractive in- 
teractions, the Feynman condition leads to a novel state. 
The two components, rather than rotate together at the 
drive frequency, rotate at angular velocities f2i,2 inversely 



FIG. 1: Schematic diagram of bound vortex pairs (for species 
with masses mi and m2 < mi) in the rest frame of the vortex 
lattice in the limit of large interspecies interaction in which 
the two vortex lattices are locked. The Magnus force Fmag 
opposes such locking and is balanced by a restoring force Frstr 
due to the interspecies vortex-vortex interaction. 



proportional to their masses, toi,2, such that: 
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while the vortex lattices of the two components lock at 
the drive angular velocity fl^, lying between Qi and Q2 
(see Fig. [ij . Qualitatively, the attractive interspecies in- 
teraction leads to an attraction between vortices of the 
two flavors. If it is sufRciently strong, vortices pair, and 
the lowest-energy vortex configuration then occurs when 
the vortex lattices of the two flavors are "locked" to- 
gether, rotate at the same rate, and have essentially the 
same density. Eq. ([T|) then reads: 

> 771if2i m2^2 
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where pl'^ are the vortex densities of the two flavors. 
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This state is strongly related to experiments in Ref. 
where a vortex lattice was locked to an optical lattice 
with a similar periodicity. 
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As we show below, this state survives within a finite 
disk about the center of the rotating condensate. The 
relative motion between the vortices and the condensate 
gives rise to a Magnus force that opposes the interspecies 
vortex attraction. Beyond a critical distance the Magnus 
force (Fig. [1]) becomes larger then the maximal pairing 
force, and the vortices become unbound. In this region, 
the two condensates and their vortices rotate together at 
the drive frequency f2^; the vortex densities in the two 
flavors are no longer equal, but reflect themselves the 
mass ratio: p\lm\ = p1/m2- Below we derive the char- 
acteristics and conditions for the vortex locking state. 

Energetics. The energy of weakly interacting Bose- 
Einstein condensates is well-described by the Gross- 
Pitacvskii functional for the condensate wave function 
-00 (a — 1, 2) for the two atomic species We consider 
the situation in which the two condensates are stirred at 
the same rate Transforming to the rotating frame, 
our problem becomes time independent. The energy of 
the two-component system in the rotating frame is given 
hy E = El + E2 + E12 , where 




+ ^n^-;if7„c(-*^)^o) (4) 

describes the energy for each atomic species and ga is the 
intraspecies coupling for bosons of flavor a. The inter- 
species interaction is given by 

£^12 =512 J dV[ni(r)n2(r)], (5) 

where Ua = IV'aP is the density of flavor a, Vtrap is the 
external trapping potential, and the z-component of the 
angular momentum operator is Lz = -~i-§^ (where (p is 
the azimuthal coordinate). 

The energy E^ of one BEG component is minimized 
via the nucleation of a vortex lattice rotating with the 
external drive f2i,. In the following, we assume that the 
coherence length £„ = \ — — — - is much shorter than 
the characteristic distance between vortices. Each vortex 
is then well-described by a small core region of size , at 
which the superfluid density drops to zero and its phase 
field accounts for the fiow around the vortex. The vortex 
core region gives rise to a small constant energy, and we 
can account for the phase field by writing i/'a = y'^e**" , 
where 9a determines a lattice of vortices with unit wind- 
ing number at the positions {rf}. We assume that we 
can write 9o. as a sum of the different vortex contributions 
9a ~ 9i + 92 + ■ ■ ■ ■ 

With these assumptions, Ea can be written in terms 



of the positions of the vortices as 

(6) 

where we have dropped terms that do not depend on the 
positions of the vortices and have also neglected effects 
due to the nonuniform superfluid density [lO|. The first 
term in Eq. ^ is the usual logarithmic interaction be- 
tween vortices, and the second is the centripetal energy, 
reflecting the fact that vortices toward the edge of the 
cloud carry less angular momentum relative to the center 
of the cloud. In a single-component rotating BEG, the 
balance of the two terms gives the Feynman condition 
([T|) . The equations describe charged particles interacting 
in two dimensions with a uniform background charge of 
opposite sign. 

The energy E12 arising from the interspecies interac- 
tion energy is less straightforward to evaluate. Unlike the 
intraspecies logarithmic interaction, this nonuniversal in- 
teraction depends on the details of the short-distance 
density variations around the vortex cores. For instance, 
in [isj to study the interaction of a vortex with an opti- 
cal lattice, a step function having the width of the BEG 
coherence length was taken. In this work we take a Gaus- 
sian depletion around the vortex core: 

n„(r) =n^(l-e-|'--''l'/«') (7) 

so that the system will be amenable to analytic treat- 
ment. For a single vortex this depletion gives the correct 
behavior at short distances, but not the long distance 
behavior, in which the density due to a single isolated 
vortex heals as ^ jr^ ■ As we show in the Appendix, the 
combined density variations due to the vortex lattice on 
scales larger than the inter-vortex separation combines 
only to change the chemical potential (which is propor- 
tional to the density correction), by ^Tr^yO^r^, which 
just reflects the kinetic energy associated with uniform 
rotation of the condensate. This can be shown to have a 
negligible effect on the vortex pairs, which are the focus 
of this work. Neglecting this piece of the density fluctu- 
ation is also consistent with neglecting the intraspecies 
core-core interactions which is a standard approximation 
0. Indeed, the short distance region is the relevant one 
for vortex locking; two-species vortex pairs become un- 
bound once their separation is comparable to the coher- 
ence lengths, where the approximate form we take for the 
density proflles is still valid. Evaluating the interaction 
integral in Eq. ([5]), and keeping only the contributions 
due to the interactions between pairs of vortices between 
different species, gives 

2 2 k"^-r^l^ 
E,2 = 3i2njng7r-|^ ■ («) 
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We note that to prevent phase separation, the criterion 
I5152I > I512P niust be satisfied. Equations ([6]) and dH) 
now give the energetics of the system purely in terms of 
the positions of the vortices. 

Forces. An understanding of the locked phases can be 
obtained by considering the forces acting on the vortices. 
Equation © leads to the well-known Magnus force acting 
on a vortex of species a 



(9) 



where k is a unit vector pointing out of the plane, Hq is 
the equilibrium superfluid density for species a (evalu- 
ated away from the vortex core), Ygp = -^VOa is the 
superfluid velocity for species a (with the vortex on which 
the force operates excluded from 9a), and is the ve- 
locity of the vortex. On the other hand, the force arising 
from the energy in Eq. ([8]) provides an attractive force 
between two vortices of different species. It has the form 
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where d is the displacement vector between vortices. 

Let us first briefly consider the unlocked case where 
the vortex interspecies interaction force is small. Since 
the force counteracting the Magnus force is too small, to 
bind vortex pairs, we must have Fmag = for any iso- 
lated vortex, which implies that the superfluid velocity 
must be the same as the vortex velocity. That is, the vor- 
tex lattice rotates with the superfiuid. Thus, because the 
two vortex lattices rotate at the same frequency, the two 
superfluids rotate together at that frequency. Accord- 
ingly, for this case, the vortex densities are not equal: 

We next consider the other extreme, in which the two 
vortex lattices are locked. Our approach is to consider 
the forces acting on a bound pair of vortices at dis- 
tance r from the center of the trap (sec Fig. [T]). As 
stated before, the Magnus force for the locked state is 
nonzero because the superfluids are rotating at differ- 
ent rates. Balancing the forces gives F^ag = -^rstr- Be- 
cause the restoring force acting on either species has the 
same magnitude, we obtain F^^^^ = ^mag- Noting that 
Va = and Vv = ^lyf (we are assuming that r is much 
larger than the distance between the two vortices) gives 
nj!j(r2„ — ill) = ^0(^2 — ^v)- This, along with the con- 
dition mifli — 1712^12 (from pi = p'^) gives the following 
relation between the angular velocities: 
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Note that unlike the restoring force, the Magnus force 
grows linearly with distance from the center of the con- 
densate. Thus, at some critical distance r,,, the pairs 
of vortices invariably become unbound from each other. 



for radii r > Vc the unlocked phase applies, and after a 
short healing region, the two condensates rotate at the 
same frequency. An expression for can be obtained by 
equating the Magnus force with the maximum possible 
value for the restoring force: 
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Note that p2|) diverges when the masses are equal. In 
addition, the bound pairs of vortices are pulled further 
apart at increasing distances from the center of the con- 
densate. The interspecies vortex separation satisfies 
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which is valid for r < Tc- This introduces a small cor- 
rection to the vortex density and creates a small shear in 
the motion of the two condensates. 

For the vortex-locked state to be stable up to the crit- 
ical radius rc, the superfluid velocities in the rotating 
frame of the vortex lattice, [v'^p — Vi,|, must not exceed 
the critical velocity of the superfluid. Otherwise, it would 
be possible to create elementary excitations from the flow 
of the superfluid around the vortices. For the two cou- 
pled superfluids, the Bogoliubov elementary excitations 
are given by 
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{n,)^ = -{El+El) ±-^{El - El)^ + lQgl,nlnlelel , 

(14) 



V(eg)2 -t- 2gan<^el and e^^ = It is 

then straightforward to compute the critical velocity Vc = 



where E"^ 



min/;(^); one obtains 
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The superfluid velocity of species 1 or 2 in the vortex 
lattice frame evaluated at (where it is maximal) is 
given by 
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^SF 



I512 



{1,2} 



/2i5{2,i} 2m{2,i} 



(16) 



The condition |vg^ — Vt,| < must be checked so that 
the vortcx-locked state is stable against the creation of 
elementary excitations. For instance, it can be shown 
that the system is stable against creating such elementary 
excitations if the conditions < ^ < 10 and 4t < 

iU — — ^ — 

i2i < 10 are satisfied. 

Numerical Simulations. Now that the expected types 
of phases have been discussed, we minimize the total en- 
ergy E ^ El + E2 + E12 as a function of the vortex 
positions given by Eqs. (O and The ability to com- 
pute analytical expressions for the gradients of the energy 
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FIG. 2: Vortex lattice for two-component condensate (with 43 vortices of each species) for different values of 7 = Igabino / (hQ^). 
Circles are shown for the theoretical prediction for the critical radius Eq. (|19p at which the vortices become unbound. 



as a function of the vortex positions V^^ijE allow us to 
apply the steepest descent method for the minimization. 
Specifically, starting with an initial configuration for the 
vortex positions {fj-"''}, we perform the one-dimensional 
minimization of 



E 



- AV{r,}-B 



(17) 



over A, where W^^ijE is evaluated at {r"}. The new 
vortex positions are given by 



{r("+^)} = {rn-A 



(18) 



and the above procedure is repeated until it converges. 

To simplify the analysis, we restrict our attention to 
the case in which the equilibrium densities and heal- 
ing lengths of the two condensate components are equal: 
no = rip = riQ, ^ = ^1 = ^2- Motivated by the example 
of a "3cs-87Rb condensate [11], we fix the mass ratio 
to be TO1/TO2 = 1.5. The vortex interaction strength is 
parametrized by 7 = ^fj^i which we vary while keep- 
ing the quantities ^^'^ for a = 1,2 fixed. (The total 
energy has been scaled by h^lyTrno^^^ .) We set the ra- 
tio of the "average" vortex lattice constant to the coher- 
ence length aiat/C==10 (which is consistent with typical 
experiments). We define aiat by = -^^^^j where 

rh — 2EL1II12. and consider a system with 43 vortices of 

mi-t-m2 ' 

each species. The results for such a calculation are shown 
in Fig. O We also plot our prediction for the critical 
radius at which the vortex pairs become unbound [see 
Eq. (fT2|) ]. which for equal densities is 



7 mi + 1712 
^\fe mi ~ 1712 



(19) 



This prediction agrees quite well with our numerical re- 
sults. 

Experimental realization. A very promising candidate 
for the realization of these locked states is a ^"^^Cs-^^Rb 



condensate mixture 'l^ , which has a mass ratio of about 
1.5. One has ex quis ite control over the self-scattering 
length of Cesium 17|, and a Cs-Rb mixture is also ex- 
pected to exhibit interspecies Feshbach resonances. This 
allows one to control the interaction gi2 over a wide 
range; such interspecies resonances have recently been 
identified for Li-Na [18| and Rb-K [ll] mixtures. 

Conclusions. In this paper, we described a novel state 
of rotating interacting condensates with unequal masses 
in the Thomas-Fermi regime. The possibility of lock- 
ing the two individual vortex lattices yields a remarkable 
demonstration of the nonintuitive behavior of superflu- 
ids: the two gasses, rather than equilibrating to the same 
speed, prefer to move at different angular velocities that 
are inversely proportional to their masses. The vortex- 
locking of the different-mass condensates is also an ex- 
ample of synchronization: a phenomenon that is ubiqui- 
tous in physics, biology, and other fields [2(|. Already 
in single-mass mixtures, a rich variety of vortex dynam- 
ics arises from the extra degrees-of-freedom, resulting in 
such effects as the formation of square vortex lattice, as 
well as topolo gica lly nontrivial defects such as skyrmions 
or hedgehogs 3, 21 1 . To investigate these effects in the 
different-mass mixtures, as well as to better establish the 
locked state we proposed in this manuscript, this system 
must be numerically investigated by solving the appropri- 
ate dynamical Gross-Pitaevskii equations. Such numeri- 
cal investigation will also allow one to find the preferred 
lattice geometry of the locked vortex-lattice. Other di- 
rections for future study involve dynamical aspects such 



as Tkachenko modes 221) of the locked state, as well as 



whether a similar state survives in the Landau regime of 
large vortex density. 
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APPENDIX 

In this appendix wc discuss the change in the supcrfluid 
density as a result of a vortex lattice, and its effect on the 
validity of approximation ([7]) and the resulting expression 
for the interspecies vortex attraction, Eq. ([8]) and ([5|). 

First consider a single species which has the energy 
functional 



E 
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d\ — I + FtrapMIV-P + Xffl^l^ . (20) 
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We write ijj — /e'^ where / is real and contains infor- 
mation about the positions of the vortices as 6* = 6*^. 
By varying / we obtain the equation determining the 
supcrfluid density which minimizes E 

|-VV + ^\'^0\\f + Vtrap/ + gf = /i/ (21) 



2m 



2m 



for the particular vortex configuration. First let us con- 
sider a single vortex taken to be at the origin VO = £ x 2^. 
The long-distance behavior r 3> ^ is obtained from 
the Thomas- Fermi approximation (neglecting the V^/ 
in Eq. BTjl and we obtain 
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where we have neglected the contribution from the trap- 
ping potential. This implies that the suppression of the 
density is due to the kinetic energy in the supercurrent, 
which counters the condensation energy of the BEG. 

Eq. (|22p seems to imply that at large distances from 
a vortex core, the interspecies vortex-vortex interaction 
will include a persistent power-law component, and die 
off only as rather than as an exponential. But the 

observation that the power law decay reflects the current- 
induced superfluid suppression, allows as to ignore the 
power-law decay in a many-vortex situation, with the 
argument as follows. Let us consider a vortex lattice, 
and evaluate / at a position which is several coherence 
lengths away from any vortex. This allows us to invoke 
the continuum approximation and write 
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where is the density of the vortices. Each individual 
contribution alone in this sum would yield a oc ^ depen- 
dence in the density corrections, but the vector sum of 
the velocities of all vortices squared is not a simple sum 



of the corrections. Inserting this into the equation 
for the density profile one finds: 
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+ + ^trap/ + gf = M/- (24) 



Thus the combined vortex effect re-normalizes the trap- 
ping potential and does not need to be explicitly taken 
into account. To get an estimate for the magnitude of 
such a renormalization, one can compare this term with 
the chemical potential 



O/at a-lat 



(25) 



where aiat is the vortex lattice constant which is small 
for typical experiments. 

For a two-component BEG, the situation is similar for 
vortices which are many coherence lengths away from 
each other. On the other hand, when the cores of the 
different types of overlap, their interaction needs to be 
explicitly calculated, and the continuum approximation 
cannot be used. This is the case for paired- vortex config- 
urations. Since the combined effect of far-away vortices 
on a locked pair is small, the locking depends only on 
the short distance density profile taken. Had we used 
the step-function potential interaction between two vor- 
tices of Ref. [l5| the results would only differ from the 



Gaussian depletion Eq.[7]by small quantitative amounts. 
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